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ABSTRACT 

We estimate the heating of the intergalactic medium due to shocks arising from struc- 
ture formation. Heating of the gas outside the collapsed regions, with small overden- 
sities (^ ^ 200) is considered here, with the aid of Zel'dovich approximation. We 
estimate the equation of state of this gas, relating the density with its temperature, 
and its evolution in time, considering the shock heating due to one-cr density peaks 
as being the most dominant. We also estimate the mass fraction of gas above a given 
temperature as a function of redshift. We find that the baryon fraction above 10^ K at 
z = is ~ 10%. We estimate the integrated Sunyaev-Zel'dovich distortion from this 
gas at present epoch to be of order 10~^. 

Key words: cosmology: theory — large-scale structure of universe — intergalactic 
medium — galaxy:formation 
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1 INTRODUCTION 

It has become evident from recent numerical simulations 
that a significant fraction of the baryons in the universe 
reside in the warm-hot phase of the intergalactic medium 
(WHIM), with temperatures of order lO'^-lO^ K (Cen & 
Ostriker 1999, hereafter C099). The gas in this phase is 
raised to high temperature by shock heating due to forma- 
tion of structure. Recent simulations by Dave et al. (2000) 
and Croft et al. (2000) have also calculated that the equation 
of state of this phase is approximately p ocT. 

There have been a few analytical attempts too to un- 
derstand the heating process in the intergalactic medium 
through analytic means. Pen (1999) pointed out that there 
is a need for non-gravitational heating in the intra-cluster 
and intergalactic medium (hereafter IGM) to avoid the con- 
straints from the soft-X-ray background. Extra heating de- 
creases the amount of clustering of the gas and therefore re- 
duces the fiux of soft X-ray radiation. Wu, Fabian & Nulsen 
(1999) have also addressed the question with detail calcula- 
tions and concluded the same. 

These are, however, relevant for heating of the gas which 
are already within collapsed objects. For example, the work 
of Wu et al (1999) refers to the heating of the gas which 
is already within a collapsed halo, with overdensities larger 
than ~ 200. The numerical simulations of C099 and Dave 
et al. (1999), on the other hand, points out the heating in 
the gas which have overdensities much smaller than this. 

It is interesting to note that Zel'dovich and his col- 
leagues had reached similar conclusions to that of the re- 



cent numerical simulations in the context of their study 
of formation of pancakes. As a byproduct of their study 
of the formation of large pancakes, they had worked out 
the magnitude of gravitational heating of the intergalactic 
gas. Although much of the earlier motivation has been lost 
now, a substantial part of their work sounds prescient. To 
quote from Sunyaev and Zeldovich (1972) — "/i is possible 
that a significant fraction of the intergalactic gas (10-50 %) 
was not subjected to compression in the 'pancakes' and was 
heated only by the damped shock waves moving away from 
them." This is exactly what the numerical simulations have 
unearthed, namely, the heating of the gas which with over- 
densities smaller than ~ 200, outside the collapsed region 
but worked upon by Shockwaves due to gravitational col- 
lapse. 

In this Letter, we attempt to understand the heating of 
this phase of IGM with the help of Zel'dovich approximation. 
Since the gas in warm-hot IGM is only mildly non-linear, this 
approximation can shed light on the gravitational heating 
process, if used within its limitations. Below, we attempt 
to estimate the amount of the gravitational heating, and 
the sate of the gas, by including other heating and cooling 
effects. We also attempt to estimate the mass fraction of 
baryons which are affected by this heating as a function of 
redshift. 

We assume a cosmological model with Qa = 0.7, fim = 
0.3 and h — 0.65, with Qb h^ = 0.015, the big bang nucle- 
osynthesis value. 
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Figure 1. The maximum temperature from slioclcs due to struc- 
ture formation is sliown as a function of redshift (solid line) . The 
dotted line is that from C099. 



2 SHOCK HEATING IN THE VICINITY OF 
COLLAPSED OBJECTS 

Consider the gas surrounding a high density peak. As the 
gas flows inwards, it is compressed, and depending on its adi- 
abatic exponent, it stops at a place away from the centre of 
accretion, and a shock wave travels outward. We will concen- 
trate on this Shockwave as it compresses and heats the very 
outer parts of the collapsed region. Sunyaev & Zel'dovich 
(1972; hereafter SZ72) tried to model this shock wave in the 
context of one-dimensional collapse of gas onto pancakes. In 
their idealized picture, as the gas flows towards the inner 
region, a singularity appears and a shock travels outwards 
through the gas. This shock velocity can be easily deter- 
mined for the perturbation of a given lengthscale A(= 2-k /k), 
assuming a single sinusoidal perturbation. Suppose the sin- 
gularity appears at a redshift Zc- They defined a parameter, 
^ which corresponds to a given Lagrangian coordinate, and 
is given by 2i2JEi£ — ^+^ , The parameter /i, therefore, is 
equivalent to a time parameter. In the case of a sinusoidal 
perturbation, it also gives the fraction of matter that has 
passed through the shock wave up to a given moment. 

The velocity of matter falling onto the shock, Vs, as 
derived by SZ72, can be generalized for any cosmological 
model as. 
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where A is the comoving lengthscale of perturbation. The 
temperature behind the shock wave is given by Ts ~ 
V]fmp/(6fcs), where rrip is the mass of a proton and fc_g is 
the Boltzmann's constant (SZ72, eqn(2)). 

Note that this temperature reaches a maximum at 
H ~ 0.5, when approximately half of the matter has passed 
through the shock. This happens when 1 + z '^ (2/7r)(l-f 2;c). 



The maximum temperature is given by, (noting that ^ 

H{z){l + z)) 

J- max 



6kB ^ ' V27ry ^^ ' \1 + Zc^ (l + 2c)2 

6kB ^ ' (l + z^r 
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where we have written L;„ ~ 1/fc for the comoving length- 
scale of the perturbation, in the notation of C099. This 
is the typical length of perturbations that goes non-linear 
at 1 + Zc- It is interesting to compare the Eqn (4) C099 
with this equation. They derived a value oi K — 0.3 from 
their simulation where the maximum temperature or, equiv- 
alently, the maximum sound velocity was given by Cf = 
KH^{Lin/{l + 2c))^. Comparing with the above expression, 
we obtain, K --^ 5/18 for a monoatomic gas. 

There is, however, a crucial difference. The parameter 
Lin in C099 is defined as the perturbation that goes non- 
linear at a given redshift z. This provided the maximum 
temperature reached by the gas at a given z. In the above 
formulation, however, there are three important epochs: zin 
is the epoch when the perturbation has an overdensity larger 
than unity and it goes non-linear, Zc is the epoch when 
the singularity appears and Zm is the epoch when fi ~ 0.5, 
when the maximum gas temperature is achieved. Naturally 
zin > Zc > Zm- Here we have assumed that zi„ --^ Zc but that 
it is larger than Zm- This difference becomes non- negligible 
especially at high redshifts. This difference is shown in 
Figure 1, where the maximum temperature reached at a 
given redshift is plotted for the A cosmological model (as 
in C099). The solid line is the prediction from the above 
formulation (with zin ^ Za> Zm and the dotted line is from 
C099. The only difference is that we have treated L;„ as 
the lengthscale of the perturbation that goes non-linear at 
a redshift slightly larger than Zra to give the maximum tem- 
perature at Zra- Naturally, there is a difference in these two 
epochs, of order of a Hubble time. 



3 EVOLUTION OF TEMPERATURE 

In order to calculate the evolution of temperature of this 
gas, we will need to take into account all sources of heating 
and cooling and their rates. Firstly, let us consider the rate 
of heating due to these shocks. We should note here that 
the gas with different initial density (or overdensity) will go 
through the Shockwave at different point of time. Gas with 
larger initial density will be closer to the singularity and 
will pass through the shock wave sooner. In the ideal case, 
the gas density profile is smooth and shock travels through 
it equally affecting all parts of it. In reality, however, we 
expect the gas closer to the collapsed region will be shocked 
and heated more effectively than the gas farther away. 

We note that the average density profile of the gas 
away from the collapsed region is expected to be of the type 
p oc r^^'"*, where r is the perpendicular distance from the 
pancake or filament (Zel'dovich 1970). To treat the differ- 
ential effectiveness of shock heating analytically, we shall 
define a timescale for the passing of the Shockwave through 
the gas as, 

A l/p^-3/2 ^3^ 



l+^c t/ 
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where p is the average ambient density, and where U — Vs/3 
is the velocity of the shock front relative to the plane of sym- 
metry (whereas Vs is the velocity of matter impinging on the 
shock (e.g., Jones, Palmer & Wyse 1981)). The extra fac- 

'pN-3/2 



tor I ^j ' then accounts for the fact that gas with larger 

overdensity is heated more effectively than that with lower 
overdensity. This prescription is valid only for collapsing gas 
with different overdensities in a given perturbation. In this 
work, we attempt to calculate the equation of state of gas in 
a given 1-a perturbation. We then estimate the equation of 
state of gas in the IGM in general using the relevant filling 
factor of such perturbations. 

We then write the shock heating rate as simply 
Tshock/t shock- For a universe with Vlf^ + fio = 1, one has, 

''^'lTT^(«A + (l + ^)'no) 
(,t ~r 2c J 



^-^ shoe 

dz 



-0.4 X 10' 



(4) 



We should, however, remember the maximum tempera- 
ture that the gas can be raised to by the shocks, as discussed 
in §2. We therefore put an upper limit on the temperature, 
as given by eqn(2). This will reflect the physical fact that al- 
though lower density gas is not shocked as effectively as the 
higher density gas closer to the filament or pancake, higher 
density gas is not heated to indefinitely higher temperatures 
this way. 

As was pointed out in SZ72, a useful approximation for 
^ is /I ~ -^(6(1— ^^^ ))''"■ We have used this approximation 
in our calculations below. 

Before discussing other sources of heating, we should 
note here that this formulation is adequate only for a limited 
duration. Although in principle, the gas infall continues until 
/i = 1, in reality, the approximations used to calculate /i 
breaks down for large values of /i. We therefore consider the 
evolution of the temperature only until /i = 0.5. 

The second heating source, which is the adiabatic com- 
pression of the gas, is easily described as, (for nb/ri;, S> 1) 



dTad 



T d{nb/nb) 



(5) 



dz 3 (nb/rib) dz 

We characterize the growth of the overdensity by the follow- 
ing equation. 



d{nb/nb) (rib/nb) 



(6) 



dz ' 1 + z 

where rj equals unity for the linear regime in a n = 1 uni- 
verse. In the quasi-linear regime, rj could be large. For ex- 
ample, the overdensity evolves as 5 oc (1 + z)^ ■^"' for the 
range of the scales where the power spectrum has n = — 2 
(Peacock 1999). We adopt a value of ?; = 2. The final result 
is found not to depend on its value strongly. 

One cooling process is due to the expansion of the uni- 
verse, and is given by, 

dT^ = JZl (7) 

dz 1 + z' ^ ' 

Cooling due to free-free radiation is given by. 



''^"=0.22nbhi^''/''^^^'''^' + '^^' 



(8) 



dz "■ (Ao + (1 + 2);^0„)V2 • 

Although cooling due to inverse Compton scattering be- 



comes important at high redshift, a,t z — 4, the cooling time 
for Compton cooling {^ 9xl0^'\l + z)-^ yr = 1.5 x 10^° yr) 
is larger than that for free- free cooling for a gas at 10^ K and 
with an overdensity of ~ 100 (~ 10^ yr). It is shown below 
that at high redshifts {zc 5; 3.5) shock heating contributes 
to the equation of state only for gas with large overdensities, 
of order ~ 100. Since Compton cooling is not as efficient as 
Compton cooling for this gas, we neglect it in our calcula- 
tion. 

Combining all the heating and cooling processes, one 
has for the evolution of the gas temperature. 



dz 



dTs 
dz 



dTai 



dT„. 



+ 



dT, 



ff 



(9) 



dz dz dz 

We present the numerical solution of this equation below. 



4 RESULTS 

The process of heating due to structure formation is essen- 
tially statistical in nature. To track it analytically, however, 
we focus on the 1-a fluctuations. These are the perturbations 
that dominate the heating at a given redshift, as found in 
C099. Fluctuations with higher degree of non-linearity at a 
given epoch would involve gas with very large overdensities, 
and does not concern us here. As has been emphasized ear- 
lier, here we are concerned with gas which is rather on the 
outskirts of highly non-linear perturbations. 

To calculate the equation of state of gas at a given 
epoch, we therefore find out the relevant lengthscale, given 
the power spectrum which is COBE normalized. For the 
physical state of the gas at a redshift z, we determine the 
perturbation which goes non- linear at Zc, so that after evolv- 
ing by a time period equivalent to /u = 0.5, we reach the 
epoch z. In other words, to find the state of the gas at 
z = 0(1, 2), we adopt Zc = {1 + z) * (7r/2) - 1 - .5(2., 3.5). 
We find the scale of the perturbations which are one-cr at 
this epoch, according to the relevant power spectrum. The 
value of Li„ at Zc ~ 0.5(2., 3.5) is ~ 6(1.6, 0.65) /i"^ Mpc. 
The initial values of Pb/p is taken to be in the range of 1-50. 
The initial temperature has been fixed at 3 x lO'* K, which is 
the temperature reached by the IGM gas through photoion- 
ization heating. The results for the state of gas at z = 0, 1, 2 
are shown in Figure 2. The dotted line shows the equation 
of state as found by Dave et al (2000) . 

We can estimate the fraction of the mass that has a 
given temperature in the following way. First, we note that, 
in the formulation of SZ72, if the density of gas which has 
just entered the shock front is pi, then the fraction of mass 
(/) that has already gone through the Shockwave is given 
as, (pi/p) ~ S/{iT^ f^). If T is the temperature of the gas 
corresponding to the initial density pi, then / is the fraction 
of mass with temperature larger than T. To be precise, this 
approximation is valid for the case of instantaneous cooling, 
which is a reasonable assumption for p ^ 0.1 {zc ^ 0.5, 
as can be seen from eqns (4) & (9)). In other words, this 
approximation is valid only for small values of f{<^ 1). 

We multiply this fraction with the (comoving) number 
density of the one-a peak, as given in Bardeen et al. (1986), 
to derive the fraction of mass which has temperature larger 
than a given value. The results for the fractions at z = 0, 1, 2 
are shown in Figure 3. Unfortunately, the limitations of the 
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Figure 2. Temperature of the gas is plotted against its overden- 
sity n^,/fn,. The solid curves are for gas at 2 = 0, 1,2 from top 
to bottom. The dotted line shows the equation of state found by 
Dave ct al. (2000) (their Fig 6). 



single sinusoidal wave approximation do not allow us to draw 
a full curve, as one has to stop at / ~ 1 (in reality, the 
above approximation is valid only for small values of /). 
The curves, however, can serve as pointers to what one-cr 
peaks can do to the diffuse IGM. 

We also show by stars the mass fraction derived by Croft 
et al. (2000) at 2: = 0, 2, and the fractions are 41% and 
5% respectively. To compare these numbers with the curves 
in Figure 3, we should remember that the curves show the 
result of heating by one-a density peaks only. In reality, 
there will be contribution from higher sigma peaks, taking 
the gas to higher temperatures at rarer places. Although 
we do not have the fraction for gas above 10^ K at z = 0, 
a naive extrapolation of the existing curve above 10® K is 
consistent with this value. We note here that the curve for 
z — Q shows that the fraction of gas above 10® K is of order 
10%. 

We note here that the mass fractions from the simula- 
tion of C099 are much larger than these. Since the simu- 
lations use different techniques and employ different resolu- 
tions, it is not obvious to what do these discrepancies owe 
their existence and if they are of much importance. 



5 DISCUSSIONS 

The equation of state for low density gas in Figure 2 de- 
pends on the particular form of the timescale for passing 
of the Shockwave in eqn (3) . as is evident from the analyt- 
ical solution (eqn (10)). In reality, modelling the efficiency 
of shock heating might add additional parameters and the 
resulting equation would therefore have some scatter. 

In Figure 3, we also plot the limits from the observa- 
tions of the soft X-ray background. After the subtraction 
of the discrete sources, it now appears that a flux of 4 keV 
cm~'^ s~^ sr~^ keV~^ at 0.25 keV can be taken as an up- 



Figure 3. The fraction of shock heated diffuse gas above a given 
temperature f2fl(> T) is plotted against logarithm of T. The 
three solid curves from right to left are for gas at 2 = 0, 1, 2 re- 
spectively. The stars arc the mass fraction at 2 = 0, and 2 = 2 
found by Croft ct al. (2000, their §2). The horizontal dotted line 
refers to the total baryonic gas in the universe and the vertical 
dotted line shows the lower limit of the temperature due to pho- 
toionization. The arrow is from the limit of soft X-ray emission 
for gas with n;,/?!;, ~ 100 at 2 = 0. 



per limit to any possible contribution from diffuse matter 
in the IGM (e.g., Wu, Fabian & Nulsen 1999). We apply 
this limit to our result for IGM at z = 0. If the gas at tem- 
perature T (with the corresponding density n{T), as given 
in Fig 2 for 2 = 0) has a filling factor er, then the flux 
at 0.25 keV is proportional to Flux {keV/cm^ s sr keV) oc 
€Tn[TfT-°-'^ exp(-0.25fceF/fei3T) {c/Hq) and this allows 
us to calculate er as a function of T. Since er is found to de- 
crease very quickly with T, we can approximate e(> T) ~ er 
and use it to put a limit on the mass fraction of gas with 
temperature larger than a given value T. This limit is shown 
as a dashed arrow in Figure 3 for gas with nb/rib ~ 100 at 
T = Tmax for 2 = 0. The limit for lower density gas is 
less restrictive. The above limit is obviously uncertain to 
the extent that the appropriate path length differs from the 
approximate c/ Ho. 

In this calculation we have used a metallicity of 0.01 so- 
lar abundance, which is the metallicity found in the Lyman- 
a absorption systems at high redshift, and is relevant for 
the diffuse gas considered here. It is possible that the metal- 
licity of this diffuse gas increases in time, as found in the 
numerical simulation by Cen & Ostriker (1999b). It is, how- 
ever, not yet known from observations how this metallicity 
evolves in time, and its dependence, if any, with the den- 
sity and dumpiness of the gas. Any metallicity will make 
the above limit more stringent; in other words, make the 
contribution of the diffuse shock-heated IGM gas towards 
to the soft X-ray background much more prominent. As the 
detail numerical simulation of the soft X-ray background 
radiation by Croft et al. (2000) shows, the contribution of 
the (enriched) shocked heated diffuse gas is comparable to 



© 0000 RAS, MNRAS 000, 000-000 



Heating of the intergalactic medium due to structure formation 5 



the upper limit of the unresolved X-ray emission in the soft 
band. 

Given the relation between density and temperature, 
we can estimate the resulting Sunyaev-Zel'dovich distortion 
on the cosmic microwave background. Defining fios to be 
the fraction of the line of sight going through hot filamen- 
tary structures, one can estimate the integrated Compton y 
parameter aas. 

The (comoving) volume fraction of 1-cr peaks from Bardeen 
et al. (1986) is of order 3%. Another way of estimating this 
is to use the fact that the gas with density <; nt occupies 
a lengthscale (lA ) {nb/nt)''^ , The Compton j/-parameter 
is approximated as, 

y^m — ^i-n-brib) —f^T, (11) 

where ot is the Thomson cross-section. For the case at 
z = 0, y for a single structure is found to be of order 
~ 4 X 10~* for Ub/fLb ~ 100. There will be, however, of 
order {{c/ Ho) / Lin) x fios '^ 30(/io3/0.03) such structures in 
one line of sight, and so the total distortion will amount to 
y~10-^ 



6 SUMMARY 

We have applied the Zel'dovich approximation to estimate 
the heating of the diffuse intergalactic medium by shocks 
associated with 1-a density peaks in structure formation at 
different redshifts. We are able to reproduce the equation of 
state of the warm-hot IGM found in recent numerical sim- 
ulations. We estimate the baryon fraction of the gas above 
10*^ K at the present epoch to be at least ~ O.lQb- The inte- 
grated Sunyaev-Zel'dovich distortion from the diffuse IGM 
filaments amounts to y ~ 10~®. 
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